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Abstract 

Optical microscopy and multi-particle tracking are used to investigate the cross-correlated diffusion of quasi two- 
dimensional colloidal particles near an oil-water interface. The behaviors of the correlated diffusion along longitudinal and 
transverse direction are asymmetric. It is shown that the characteristic length for longitudinal and transverse correlated 
diffusion are particle diameter d and the distance z from particle center to the interface, respectively, for large particle 
separation r. The longitudinal and transverse correlated diffusion coefficient Dj(r) and D±(r) are independent of the 
colloidal area fraction n when «<0.3, which indicates that the hydrodynamic interactions(Hls) among the particles are 
dominated by His through the surrounding fluid for small n. For high area fraction «>0.4, the power law exponent for the 
spatial decay of D±(r) begins to decrease, which suggests the His are more contributed from the 2D particle monolayer self 
for large n. 
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Introduction 

The dynamic behavior of confined colloidal suspensions has 
recently received a considerable amount of attentions [1—6]. 
Microfluidic devices, porous media, fluid interface or cell 
membrane [7-15] are real-world circumstances in which particles 
are usually spatially confined. Colloid behavior in spatially 
confined environments is more complicated than that in 
unbounded three-dimensional (3D) fluid bulk. The correlated 
diffusions coefficient D\\ and Z>i between a pair of particles in 3D 
bulk are well known as D\\,D± oc 1/r and D\\=2D± [16]. The 
hydrodynamic interactions (His) between a pair of particles 
confined by solid walls decay with the inter-particle separation r 
as 1/r 2 [4,17,18]. The behaviors of the longitudinal and transverse 
correlated diffusion of particles confined by a two-dimensional 
(2D) interface or viscous membrane are different [2,16-21]: 
D\\cc\/r and D±azl/r 2 , which is caused by the break in the 
spatial symmetry. 

The diffusion behaviors of particles reflect the influence of the 
boundary conditions [22-26]. Solid wall condition with a non-slip 
boundary [27], which could cut off the fluid field in its vicinity, has 
been studied extensively. However, there are comparatively few 
experimental studies on the influences of the fluid-fluid interface, 
which could partially transforms the surrounding flow field, on the 
correlated diffusion of particles. Ref. [28] investigated the 
influence of a fluid-fluid interface on correlated diffusion and 
presented a scaling method to yield a master curve capturing a 
transition from 3D bulk to the interface dominated. While the 
colloidal monolayer in 3D bulk approaches to the interface, the 
dynamic responses of the correlated diffusion along the two 



normal direction are asymmetric, which need to be further 
investigated in detail. 

In this paper, we report an experimental measurement of the 
cross-correlated diffusion of colloidal particles near a water- 
decahydronaphthalene (decalin) interface. It is shown that for 
larger particle separation, along the line connecting the centers of 
the particles and the direction perpendicular to this line, the 
characteristic length is particle diameter d and the distance 
between particle monolayer and the interface z, respectively. The 
influence of a fluid-fluid interface on the dynamic behavior of the 
colloidal monolayer tends to saturate, when the particle separation 
is much larger than the distance between the interface and the 
monolayer. These phenomena reveal that the hydrodynamics 
interactions (His) of monolayer particles occur through the 
colloidal monolayer self and the surrounding fluid, and the 
weights of His through these different paths (the colloidal 
monolayer or its surrounding fluid) change with the particle 
separation. In addition, we measured the correlated diffusion of 
particles for high colloidal density. In the case of colloidal area 
fraction n>0.4, D« and D± begin to depend on n, which suggests 
His through colloidal monolayer dominate even for large particles 
separation. 

Materials and Methods 

Two types of colloidal particles, both with a diameter 
d = 2.0[im, were used: polystyrene (PS) latex spheres and silica 
spheres. The silica and PS spheres were purchased from Bangs 
and Invitrogen, respectively. PS microspheres with sulfate 
functional groups on the surface are negatively charged. The 
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surface functional groups of silica sphere is Si-OH. The PS and 
silica samples were suspended in deionized water of 18.2 MQ-cm. 
By centrifugation, we cleaned the samples seven times to eliminate 
any possible surfactant before use. 

The experimental setup is shown as in Fig. la. The sample cell 
is composed of a stainless steel disk with a hole of the diameter 8.3 
mm, and the bottom of this hole is sealed with a 0. 1 -wim-thick glass 
cover slip, that also serves as an optical window. In addition, the 
inner container of the cell is divided into two layers, as shown in 
Fig. la. The cleaned water-sphere suspension was filled into the 
bottom layer, and then, decalin (a mixture of cis and trans with a 
density of 0.896 g/crrr) purchased from Sigma-Aldrich was added 
to the top of the suspension, completely filling the inner container. 
We cover the top of the inner container with another cover slip. 
We let the cell settle upside-down for two hours after flipping 
quickly. The fluids in the cell are completely pinned by the 
surrounding solid boundaries. Thus, the water layer could remain 
on the top and the oil will stay on the bottom. With the aid of 
gravity, colloidal particles fall onto the water-decalin interface and 
form a monolayer, as shown in Fig. la. All DI water was used 
within 30 min after preparation. The colloidal particles preferen- 
tially remain in the water rather than sink into decalin for they are 
hydrophilic. Particles remain in the medium (water) with the 
higher dielectric constant and the image charge in decalin has 
same sign as the surface charge of the particles [29]. Due to 
Coulomb repulsion, particles keep a certain distance away from 
the interface. The influence from cover slips can be neglected 
because the distance between the cover slip and the colloidal 
monolayer is about 0.8 mm and sufficiently large. 

Using an inverted Olympus FX71 microscope, the motion of the 
colloidal particles was recorded by a digital camera (Prosilica 
GE1050) at a rate of 14 frames per second. Each image sequence 
comprises 500 consecutive frames (taken over ~ 35.7 sec) and each 
image is 176x 176 jim 2 . From these image sequences, we obtain 
the particle positions and construct their trajectories with a 
homemade software. The center positions of particle were 
calculated using the geometric average of pixels with gray-degree 
weight. Generally, the positions could be also obtained using 2D 
surface peak fitting with a higher spatial resolution than using the 
geometric average of pixels, up to within 0. 1 pixel. Our calculation 
shows the results obtained using two methods agree within 5 per 
cent deviation. This is because the particles are highly symmetrical 
and include enough pixels (more than 100 pixels per particle). The 
real length of a pixel is approximately 170 nm. The spatial 
resolution was estimated to be third or half pixel, approximately 
60-100 nm. 

Results and Discussion 

Positions of the particle monolayers 

The PS monolayer is located at a higher position above the oil- 
water interface than the silica monolayer because the PS spheres 
are less dense and have a greater surface charge than the silica 
particles. The average distance z between the particle monolayer 
and the interface can be estimated by measuring the single 
diffusion coefficient of the particles. The single-particle mean 
square displacement (MSD) <(Av(t)) 2 > = (\s(t+x)—s(t)\ 2 } was 
calculated from the particle trajectories s(t). The MSD curves with 
typical different n are plotted in Fig. 2. According to the equation 
((As(r)) 2 y = 4D$ * t, we obtained the long time self-diffusion 
coefficient Ds for different particle area fractions n. Figure 3 
exhibits the relationship of the diffusion coefficient Ds/ Do to n for 
the PS and silica spheres, where Do =kgT '/3nt] w d is the diffusion 
coefficient of a single particle in water. We fit the results with the 
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Figure 1. (color online) (a) A schematic of the process of the 
sample preparation, (b) An optical microscope image of the silica 
spheres (d = 2.0 fori) suspended near the water-decalin interface at an 
area fraction of » = 0.01 (note that the water is on the top and the 
decalin is on the bottom). 
doi:10.1 371/journal.pone.00851 73.g001 

second-order polynomial, 

D s /D 0 = x(\-pn-yn 2 l (1) 

which is shown in Fig. 3 as the solid lines. Where cc, /? and y are the 
fitting parameters, whose fitted values are given in Table I. 

The parameter a represents the local viscosity experienced by a 
single sphere at the dilute limit, whose value can be used to 
estimate the distance between the interface and the particle 
monolayer z. The distance from a particle's center to the interface 
z can be written as [26,30]: 

z = 3(2^,-3^) 

a 16( )? „, + )7o )(a-l)' 1 ' 

where a is the particle radius, rj w and r\ 0 are the viscosity of water 
and decalin, respectively. Substituting the fitted values of a, into 
Eq. 2, we obtain the distance z = 2.3 + 0.2 fim for PS spheres and 
1.4 + 0.1 fjm for silica spheres. 

Cross-correlated diffusion coefficient D\\ and D± 

Tracking individual particle's trajectory, we obtain particle's 
cross-correlated motion via the ensemble of averaged tensor 
products of the particle's displacements [16]: 
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Figure 2. (color online) The curves of single-particle mean 
square displacement (MSD) for PS spheres and silica spheres. 

The area fractions n are 0.05 (square black symbols) for PS spheres, 0.05 
(dot red symbolsjand 0.52 (triangle green symbols) for silica spheres. 
doi:1 0.1 371 /journal.pone.00851 73.g002 



M xy (r,x) = < A^ x {t,t)A^{t,%) > w> „ (3) 

where As y (t,x) = s l y (t-\-T) — s l ^t). The superscript i and j indicate 
distinct particle, and the subscript x and y represent different 
coordinates, and r is the separation between particles i and j (as 
shown in the inset of Fig. 4). We focus only on the diagonal 
elements of the tensor product, M rr and Mgg for the off-diagonal 
elements are uncorrelated. The diagonal element M rr indicates the 
correlated motion along the line connecting the centers of the 
particles (longitudinal), and Mgg represents the correlated motion 
perpendicular to this line (transverse). It is shown that for a small 
time lag T, the measured correlated motions M„ and Mgg are 
linear functions of x. And, the cross-correlated diffusion coeffi- 
cients are defined as Z)|| =M rr /2i and Z)j_ = Mgg/2x. 

Typically, 10 5 particle pairs were used to calculate the diffusion 
coefficients D\ and D± in Fig. 3. For example, with the area 
fraction (n = 0.01), 10 short movies were acquired and each movies 
containing 500 frames. Nx90 particles were found in each frame. 
The number of particle pairs used for the averaged diffusion curve 
is [N(N-\)/2] * 10 * (500— l)/2 = 10 7 . Because each diffusion 
curve usually contains 100 data points, each data point is thus 
averaged over 10 s particle pairs. These numbers vary slightly 
among the colloidal samples with different concentrations. As 
Savin and Doyle [31] suggested, averaging of 10 4 data points is 
enough to reduce the statics error in most cases. Because the 
exposition time to capture a frame was 30 ms, and the diffusivity 
of particles during that period is approximately 50 nm. The 
dynamics error could be also safely neglected. 

Figure 4 exhibits D\\ and D± of the PS and silica spheres as 
functions of r for different area fraction n. Each curve of a fixed n 
was obtained by averaging at least 10 6 particle positions. The 
deviation of each n was given in the legend of Fig. 4. The 
behaviors of D\ and D±_ of the PS and silica monolayer are 
qualitatively similar. With the increase of r, Z)j and D± of the PS 
and silica spheres decrease in the form of power law \jr'\ In the 
case of area fraction n<0.3, the curves of D\(?) or D±(r) for 
different n almost lie on a single curve respectively, i.e., the particle 
area fraction has little effect on their cross-correlated motion. 
While n > 0.4, the power law exponent (/., called the decay rate in 




n 

Figure 3. (color online) The measured self-diffusion coefficient 
As scaled by A, as a function of the particle area fraction n for 
silica and PS spheres. Different symbols represent data for different 
particles. The solid lines show the second-order polynomial fitting 

DslD 0 = 3.(\-ftn-yn 2 ). 
doi:10.1371/journal.pone.0085173.g003 

the later part of the paper) of D\\{r) and D±(r) decrease with the 
increase of n, as shown in Fig. 4(b). 

The mechanism behind the independence of n mentioned 
above is in line with the theoretical description given by 
Oppenheimer et al. in Ref. [21]. In the particle monolayer, the 
far-field response of the correlated particle motion mainly arises 
from the momentum diffusion through the 3D surrounding fluid. 
The viscosity of 3D fluid suspension is usually a function of volume 
fraction (j) of particles. In our system (f> can be regarded as almost 
zero, regardless any change in the area fraction n of the 
monolayer. Thus, the viscosity of the 3D fluid almost keeps as a 
nearly constant. The far-field 3D His between the particles are 
independent of the area fraction n, and as a result, the change of 
area fraction n does not affect the cross-correlated diffusion D^(r) 
and D±(r). Indeed, Cui et al. [17] also observed similar 
concentration independent effects for particles confined between 
two plates, but the reasons for their phenomena are different. In a 
confined particle monolayer, the fluid momentum is absorbed by 
the solid boundaries, and hence, the far-field fluid response arises 
solely from the mass-dipole perturbations, which are not 
influenced by the presence of neighboring particles. 

Influence of the interface 

To focus on the influence of the oil-water interface on the 
correlated diffusion coefficients, first, we averaged the curves of 
«<0.3 shown in Fig.4. Then, the averaged curves D\\{f) and 
D±(r) are normalized by the single diffusion constant ocDq. We 
found that the amplitude of D\\{r)/<xDQ of the silica spheres is 
almost 1.2 times larger than that of the PS spheres, which number 
is just equal to the ratio of ft (0.78/0.62=1.2, as displayed in 

Table 1. The distance z from particle's center to the interface 
and the fitted values of the a., \i and y in Fig. 3. Here a= Ifim is 
the radius of the particles. 



Sample 


z/a 






a 






n 








PS 


2.3 


+ 


0.2 


0.87 


+ 


0.01 


0.62 ± 


0.12 


1.54± 


0.61 


Silica 


1.4 


± 


0.1 


0.77 


± 


0.01 


0.78 + 


0.02 


0.78 + 


0.03 



doi:1 0.1 371/journal.pone.00851 73.t001 



PLOS ONE | www.plosone.org 



3 



January 2014 | Volume 9 | Issue 1 | e85173 



Correlated Diffusion of Colloidal Monolayer 




10 L 



£ 10- 



Q 

d 



10"' 



10'' 



! (b) Silica 




: 








: s 8$. 






_ ■ n=0.012±0.003 






• n=0.032±0.004 






a n=0.056±0.003 


O * 




. t n=0.078±0 004 






♦ n=0.100±0.004 






•< n=0.180±0.005 






r n=0.31 1+0.007 
• n=0.403±0.007 




<] .33 OSO 
o □ 0& : 


* n=0.516±0.014 







10 u 



10 1 

r (urn) 



10' 



Figure 4. (color online) The measured correlated diffusion 
coefficient D (solid symbols) and D_ (open symbols) as a 
function of the inter-particle distance r with various values of 
the area fraction n for (a) PS and (b) silica spheres. The solid lines 
with slopes — 1 are visual guides. The geometry of measuring the 
correlated diffusion is depicted in the inset. 
doi:10.1371/journal.pone.0085173.g004 

Table I). This result is convincing because the parameter /? reflects 
the strength of the two-body His between two spheres, which 
shows that the results agree well for the measurements obtained by 
one or two particle diffusion methods. After the particle separation 
r is normalized by particle's diameter d, D^(r/d)/(fio(.Do) in 
different colloidal monolayer overlaps for large r, shown as in 
Fig. 5(a). 

In the perpendicular direction, D±(r/d)/(/iciDo) for PS and 
silica spheres deviates seriously. However, if we normalize the 
particle separation r with the distance z, D j_(r / z) j '(/taDo) of 
different colloidal monolayers overlap for large separation 
(r/z>6), as shown in Fig. 5(b). Different methods to normalize 
the particle separation r indicate that the characteristic length of 
longitudinal and transverse correlated diffusion is d and z, 
respectively, for far-field response. 

We compare our results with previous studies for the following 
boundary conditions: colloidal spheres dispersed in an unbounded 
3D bulk, confined by a solid wall or at an air-water interface [1 7— 
20] . We find that the cross correlated motion of colloidal spheres 
near an oil-water interface is similar to that at a viscous membrane 
[19,20]: D\\(r)\/r and D±(r)\/r 2 when r is large enough. For a 
viscous membrane, r called "large enough" corresponds to ry>L s ; 
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Figure 5. (Color online) (a) The mean correlated diffusion 
coefficient D||/(j3aDo) (solid symbols) and £>i/(/?xD 0 ) (open 
symbols) as a function of r/d for PS and silica particles, (b) 

The mean correlated diffusion coefficient D|/(/JaZ)o) (solid symbols) 
and D±/(J3aDo) (open symbols) as a function of r/z for PS and silica 
particles. The square represents data for the PS particles and the dots 
represent the silica particles. The solid lines with slopes — 1 and —2 are 
visual guides. 

doi:10.1371/journal.pone.0085173.g005 

note that L s = r\ m /2r} h is the Saffman length [32], r\ m is the 2D 
membrane viscosity and rj b is the 3D bulk viscosity of the 
surrounding fluid. For our system, r called "large enough" 
corresponds to r»z or ry>d in different directions. 

At first glance, the result that the decay rate of -Di(f) for the 
silica spheres is larger than that for the PS spheres seems to conflict 
with the fact that a larger decay rate should correspond to a 
smaller monolayer viscosity. Because the silica monolayer are 
closer to the interface, the silica spheres are located in a more 
viscous environment in comparison to the PS spheres. Thus, the 
decay rate of D±(r) for the silica spheres should be less than that of 
the PS spheres. In fact, the viscosity )j(z) experienced by a sphere is 
the 3D viscosity of the local surrounding fluid, not the 2D viscosity 
\] m of the membrane. The influence of this surrounding viscosity 
r\(z) has been removed by 1/ (otDrj) scaling, where a represents the 
effect of the local viscosity >?(z). The finding that the value of the 
decay rate of D^{r) for the silica is larger than that of the PS 
spheres stems from the His that are modified by the interface. 
Because oil is more viscous than water, the flow field induced by 
the particles in water is suppressed by the oil-water interface. In 
addition, the His between particles at a distance r from each other 
decay rapidly when the monolayer is closer to the interface, which 
leads to a larger measured decay rate of D±(r) for a smaller z. 
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His through three paths 

Similar to the mechanism of a membrane near a solid wall [19], 
the His among the particles in our system are transmitted through 
three paths: I) the 2D flow through the monolayer, II) the flow 
through the fluid layer sandwiched between the monolayer and 
the interface, III) the flows through the upper water layer and the 
below oil layer. Whereas the His through path I are dependent on 
2D monolayer' viscosity r\ m , the His through paths II and III are 
dependent on the 3D viscosities rj 0 and t] w of the surrounding fluid. 

The His through path I entail that the cross-correlated motion 
is a function of n; indeed, rj m is a function of n usually. The His 
through path I only contribute to the correlated motion when the 
separation r is less than the order of Saffman length 
L s = r i m /( r io + r lw) [32]- In our system, the viscosity r\ m of the 
monolayer is very small, and we do not measured the n effects for 
n<0.3. As r increases, the weight of the His through path I 
quickly decreases and the His through path II contribute to the 
correlated motion, rendering the correlated motion a function of z 
(Fig. 5). For a large separation r»z, the His through path III are 
dominant. Because rj 0 and r\ w shows little change with the area 
fraction n, the correlated motion should be independent of n (Fig. 4 
for «<0.3). In addition, the decay rate of the correlated motion of 
the particles is independent of z for r»z. In the intermediate range 
of r, the relative weights of the His through paths II and III change 
gradually with r. The decay rate of the cross-correlated motion 
shows a crossover tendency. 

For a very high area fraction w, the averaged inter-particle 
separation permits an strong interaction between the particles, 
which means that the 2D monolayer viscosity rj m would become so 
large that the order of L s would be extended to the entire 
measuring range of r [28]. Then, the cross-correlated motion 
should be a function of «, which may be supported by Fig. 4b: as 
the curves of D±(r) for « = 0.4 or 0.52 do shift upwards in 
comparison to those for «<0.3. 
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